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ABSTRACT 

Pulsar radio emission is modelled as a sum of two completely po- 
larized non-orthogonal modes with the randomly varying Stokes pa- 
rameters and intensity ratio. The modes are the result of polariza- 
tion evolution of the original natural waves in the hot magnetized 
weakly inhomogeneous plasma of pulsar magnetosphere. In the course 
of wave mode coupling, the linearly polarized natural waves acquire 
purely orthogonal elliptical polarizations. Further on, as the waves 
pass through the cyclotron resonance, they become non-orthogonal. 
The pulse-to-pulsc fluctuations of the final polarization characteris- 
tics and the intensity ratio of the modes are attributed to the tem- 
poral fluctuations in the plasma flow. 

The model suggested allows to reproduce the basic features of 
the one-dimensional distributions of the individual-pulse polarization 
characteristics. Besides that, propagation origin of pulsar polarization 
implies a certain correlation between the mode ellipticity and posi- 
tion angle. On a qualitative level, for different sets of parameters, 
the expected correlations appear compatible with the observed ones. 
Further theoretical studies are necessary to establish the quantita- 
tive correspondence of the model to the observational results and to 
develop a technique of diagnostics of pulsar plasma on this basis. 

Key words: waves — plasmas — polarization — pulsars: general 



1 INTRODUCTION 

Pulsar magnetosphere contains the ultrarelativistic 
electron-positron plasma which streams along the 
open magnetic lines. The radio emission observed 
from pulsars is undoubtedly associated with the 
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processes in the plasma flow, and its characteristics 
are believed to reflect the properties of the ambi- 
ent plasma. Whatever the radio emission mecha- 
nism, it should give rise to the natural waves of the 
plasma. The ultrarelativistic strongly magnetized 
plasma of pulsars allows two types of non-damping 
natural waves, the ordinary and extraordinary 
ones, which are linearly polarized in orthogonal di- 
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rections, in the plane of the ambient magnetic field 
and perpendicularly to this plane, respectively. 
Pulsar beam is believed to be an incoherent mix- 
ture of the two polarization modes. Since the elec- 
tric vector of each mode is related to the magnetic 
field direction, the position angle (PA) of linear po- 
larization should vary smoothly across the pulse as 
the pulsar beam rotates with res pect to an observer 
iRadhakrishnan fc Cook^ll969^ . The expected S- 
shaped swing of PA is indeed observed in a number 
of pulsars. In other cases, abrupt nearly orthogonal 
jumps of PA can break the smooth swing, testifying 
to the presence of both polarization modes. In gen- 
eral, these modes can be markedly non-orthogonal 
(e.g. 'Gil. Snakowski fc S tinebrindll99ll: iGil et all 
[W92; Gil & Lvnc 1995). 

An extensive analysis of observational data 
has led to the conclusion that the two po- 
larization modes are simultaneously present in 
pulsa r radiation jMcKirmon & Stincbring 199^, 
I2OOCII . From the theoretical point of view, this is 
the only way to account for partial depolariza- 
tion of pulsar radiation, since the natural waves 
are completely polarized by definition. Thus, the 
model of superposed polarization modes with ran- 
domly varying intensities has weighty observa- 
tional and theoretical grounds. At the same time, 
it faces serious difficulties and seems too sim- 
plified to explain diverse and complicated be- 
haviour of the individual pulse polarization. Gen- 
erally, the observed modes have elliptical rather 
than linear polarization and can be non-ortho gonal 
llKarastergiou et al.l I2OOII: iMcKinnonI l2003l) . Be- 
sides that, the observed pulse-to-pulse polariza- 
tion fluctuations cannot be solely attributed to the 
variations of mode intensities - the polarization 
characte ristics of the modes should fluctua te as 
well (e.g. lKarastergiou. Johnston fc Krameii 2003). 
These difficulties can be avoided if one take into 
account propagation effects during wave passage 
through the flow of pulsar plasma. Then the fluc- 
tuations in the observed polarization characteris- 
tics can be attributed to the fluctuations in the 
parameters of the plasma. It should be noted that 
the propagation origin of pulsar polarization im- 
plies a certain correlation between the final polar- 
ization characteristics, since they are determined 
by the instanteneous state of the plasma in the 
regions of significant polarization evolution inside 
the pulsar magntosphere. The indications of such 



a correlation are indeed discovered in obser vations 
jEdwards fc Stapperjiiool lEdwardjliooi) . 

Deep inside the magnetosphere, the plasma 
number density is so large that the natural waves 
propagate in the regime of geometrical optics, with 
the polarization planes being adjusted to the lo- 
cal magnetic fleld direction. Further along the tra- 
jectory, as the plasma density decreases consid- 
erably, the wave mode coupling starts: Each of 
the incident waves becomes a coherent mixture 
of the two natural waves peculiar to the ambi- 
ent plasma. As a result, the waves acquire some 
circular polarization and a shift in PA, so that 
it no longer reflects the magnetic field geometry. 
The analytical and numerical tracings of the po- 
larization evolution because of wave mode cou- 
pling have proved that at the conditions relevant 
to the magnetosphere of a pulsar this process 
can be efficient enough to account for the typi- 
cally observed circular polarization and PA shift 
(iLvubarskii fc Petro vall999HPetrova fc Lvubarskiil 
2000; Petrova 2003). For given plasma parame- 
ters, both types of waves acquire the same shift 
in PA and the same fractional circular polarization 
of opposite senses. Thus, the outgoing waves have 
the elliptical polarization purely orthogonal at the 
Poincare sphere. This representation is compati- 
ble with the bulk of observational data. However, 
sometimes a clear evidence is met that the polar- 
izatio n modes are non-orthogonal (e.g. iMcKinnonI 
l2003h . 

The non-orthogonality of polarization modes 
can be attributed to the cyclotron absorption 
in pulsar magnetosphere, i f this p rocess follows 
the wave mode coupling JPetrovai ■2005. . .200Q . 
A similar idea has recently been discussed in 
jMelrose et alJl2006ll . Typically the region of cy- 
clotron resonance lies beyond the region of efficient 
wave mode coupling, but they can be quite close to 
each other, in which case the plasma number den- 
sity is still large enough for the resonance to affect 
wave polarization considerably. For the two types 
of natural waves, the rate of cyclotron absorption is 
slightly different. Since the waves entering the res- 
onance region present already a coherent mixture 
of the two natural waves and these constituents are 
absorbed not identically, the polarization state of 
the waves changes. The original ordinary and ex- 
traordinary waves suffer different polarization evo- 
lution and become non-orthogonal. 
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The aim of the present paper is to study statis- 
tics of the propagation-induced polarization of pul- 
sars resulting from the fluctuations in the plasma 
flow. Special attention is paid to studying the cor- 
relation between the mode ellipticity and PA be- 
cause of propagation effects. Section 2 contains the 
general equations of polarization evolution in pul- 
sar plasma and an example of numerical tracings 
of the ellipticity and PA. The statistics of the final 
polarization parameters in the plasma with ran- 
domly varying parameters are examined in Sect. 3. 
The numerically simulated histograms of PA and 
ellipticity, as well as the two-dimensional scatter 
plots of the Stokes parameters are given. The re- 
sults of the paper are discussed and summarized in 
Sect. 4. 



2 GENERAL THEORY OF 

POLARIZATION EVOLUTION IN 
PULSAR PLASMA 

Let the radio waves propagate through the ultrarel- 
ativistic electron-positron plasma, which streams 
along the open field lines of the dipolar mag- 
netic field of a pulsar. Because of continuity of 
the plasma flow within the open field line tube, 
the plasma number density decreases with distance 
from the neutron star just as the magnetic field 
strength, N oc B oc z~^. Thus, the radio waves 
propagate in the weakly inhomogeneous medium, 
with the scale length for change in the parameters 
much larger than the wavelength. The plasma is 
assumed to be hot, so that the region of cyclotron 
resonance is also much larger than the wavelength. 

The waves are believed to originate deep in- 
side the magnetosphere, where the plasma number 
density is large enough to provide the conditions of 
geometrical optics and the magnetic field is strong 
enough for the electron gyrofrequency to be much 
larger than the radio frequency in the plasma rest 
frame. As both A'^ and B decrease rapidly along 
the trajectory, both these conditions are ultimately 
broken. First of all, the scale length for beats be- 
tween the natural waves, Lh ~ — i— , where An is 
the difference in their refractive indices, becomes 
comparable to the scale length for change in the 
medium parameters, Lp ~ z, so that geometrical 
optics is violated and wave mode coupling holds. 
In case of the cold ultrarelativistic strongly mag- 
netized plasma. 



An = 



2iji 



where ujr> 



(1) 



is the plasma frequency, 7 



the plasma Lorentz-factor, 9 the wavevector tilt to 
the ambient magnetic field and it is taken into ac- 
count that the wave propagation is generally quasi- 
transverse, 9 ^ 1/7. Then the characteristic radius 
of wave mode coupling is determined by the rela- 
tion 



(2) 



2ujl{rp)rp ^ ^ 
')^cu>9'^(rp) 

fn case of hot plasma, equations (1) and (2) remain 
the same, but 7 means some characteristic Lorentz- 
factor of the plasma particle distribution. For the 
conditions relevant to pulsar magnetosphere, the 
radius of wave mode coupling, r^, c an be estimated 
as (equation (16) in|Petrw3(20o3): 

3/2^ -3/2 1/2 nl/2 -l/2a-l 



^ = QA^p-'/^-y^r^' Bl{lu-''^9-_\ . (3) 

Ifere is the light cylinder radius, P the pulsar 
period, k the plasma multiplicity, K2 = -B* 
the magnetic field strength at the surface of the 
neutron star, i3*^2 = j^tq, ^ the radio frequency, 

that Tp is determined by the basic pulsar parame- 
ters, P and B*, as well as by the parameters of the 
plasma flow, k, and 7. Besides that, Tp strongly de- 
pends on the angle 6, which introduces a significant 
uncertainty in the estimate (3). The value of 9 may 
differ within an order of magnitude (0 ~ 0.1 — 1) for 
different pulsars and for the rays observed at dif- 
ferent longitudes in a given pulsar. On the whole, 
one can conclude that in the majority of pulsars for 
all the rays forming the pulse profile Tp lies inside 
the light cylinder. 

It should be noted that the wave mode cou- 
pling holds only if the wavevector goes out of the 
initial plane of magnetic lines. In pulsar case, this 
condition is certainly satisfied, e.g. because of the 
magnetosphere rotation. Let us choose the Carte- 
sian coordinate system with the z-axis along the 
wavevector and the x-axis in the initial plane of 
magnetic lines (see Fig. 1). The ray is emitted along 
a field line of the dipolar magnetic field at an angle 
1/; to the magnetic axis. In the non-rotating mag- 
netosphere, the ray would propagate in the plane 
of magnetic lines and at distances much larger 
than the emission altitude would make the angle 
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~ tp with the ambient magnetic field because of 
divergence of the magnetic hnes, so that bx ^ tp, 
by — 0. Because of the magnetosphere rotation, by 
increases with distance oa z/vl (the exact formulas 
for bx and by, which allow for the rotatio nal aber- 
ration , are given by the equation (3) in IPetroval 
(|2003)). The quantity /j, = {by/bx)r^ ~ rp/{rLip) is 
the key parameter which determines the efficiency 
of wave mode coupling. It should be noted that al- 
though as a rule rp/vL ^ 1, can be of order unity 
and even larger. The process of wave mode coupling 
is most efficient ai fi — 1, in which case the resul- 
tant degree of circular polarizatio n of the waves is 
as large as about 40% (cf. Fig. 5 in lPetroval j200lD ). 
Generally speaking, fj, depends on the geometry of 
ray propagation and can differ substantially for the 
rays observed at different pulse longitudes and fre- 
quencies as well as for different pulsars. Moreover, 
Vp and n may vary from pulse to pulse because of 
fluctuations in the plasma parameters (cf. equation 
(3)). 

The process of wave mode coupling holds as 
long as Li, ~ Lp. As the plasma number density 
decreases further, this process ceases and the waves 
propagate as in vacuum, preserving their polariza- 
tion state. For this reason, the wave mode coupling 
is usually called polarization-limiting effect. How- 
ever, actually it is the final stage of polarization 
evolution only if the region of cyclotron resonance 
lies infinitely far. Otherwise the cyclotron absorp- 
tion also affects the wave polarization. For the par- 
ticles with the characteristic Lorentz-factor, the 
condition of cyclotron resonance, uj~/9^ /2 — luh = 
— , is met at 

^ = o.55p-'Bli^u^'^S;lJ'ezl/' (4) 

(e.g. lPetrovalbOOeri . Similarly to the radius of wave 
mode coupling, the resonance radius is strongly af- 
fected by the uncertainty in the angle 6, but gen- 
erally Vc lies within the light cylinder. Comparing 
the equations (3) and (4), one can conclude that 
for a wide range of pulsar parameters rc ~ rp 
and the process of cyclotron absorption can play 
a marked role in determining the final polarization 
of radio waves. The quantity rj = rp/vc is another 
key parameter of the polarization evolution in pul- 
sar magnetosphere. 

Generally speaking, the plasma of pulsar mag- 
netosphere is gyrotropic, i.e. the distribution func- 
tions of electrons and positrons are somewhat dif- 



ferent. Correspondingly, the natural waves are lin- 
early polarized only in the vicinity of the emis- 
sion region, where the approximation of infinitely 
strong magnetic field is still valid. In the outer mag- 
netosphere, in the regions of significant polariza- 
tion evolution (at r ~ rp,rc) the natural waves 
of the medium are generally elliptical. Unfortu- 
nately, the question on the net current density 
and its distribution in the magnetosphere is still 
open. Given that the electrons and positrons differ 
only in the number densities and this difference is 
of the order of the Goldreich- Julian number den- 
sity, i.e. AN/N ~ the relative contribution 
of the plasma gyrotropy to the polarization evolu- 
tion in the region of wave mode coupling can be 
estimated as 0.1k,2 ^ S -I'yi.s'^' /<jJh jPetrova 2006) : 
here w' is the wave frequency in the plasma frame, 
uj' = uj^9'^/2. One can see that if the wave mode 
coupling holds close enough to the region of cy- 
clotron resonance,i.e. if ui' /luh ~ 1, the contribu- 
tion of the intrinsic ellipticity of the natural waves 
can be significant, at least for the rays of a spe- 
cific geometry. However, as the true form of the 
net charge density is unknown, we leave out the 
plasma gyrotropy, keeping in mind that this effect 
can also be significant and may somewhat alter our 
quantitative results. In the present paper, we con- 
centrate on studying the evolution of linearly polar- 
ized natural waves as a result of propagation effects 
in pulsar plasma. In our consideration, the wave 
ellipticity arises and changes purely on account of 
wave mode coupling and cyclotron absorption. An 
opposite approach has recently been developed by 
Luo & Melrose ( 2004) and Melrose & Luo ( 2004a), 
who investigated the characteristics of the ellipti- 
cally polarized natural waves at r ~ rp, ignoring 
the wave mode coupling in this region. The effect of 
cyclotron absorption on the elliptical natural waves 
has been considered in Melrose & Luo { 2004^. 

An exact treatment of polarization transfer in 
pulsar plasma with account for the wave mode cou- 
pling and cyc lotron absorption can be found in 
iPetroval ll2006l) . The evolution of the Stokes pa- 
rameters of the natural waves is described by the 
following set of equations: 

d/ 2 ^ / 

— = -77 swF2q 1 5-^ - 2^T]sF2U , 

dw \ rj^w^ J 

^ = ^r,^swF2l (1 - - 2fi7^sFiv , 

aw y r]-^w^ J 
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- — — ri swriv 
dw 



- — — —ri swriu i 
aw ' 



2iir,sF2l, 



2firisFiq . 



(5) 



Here w = r^jz 



v.p. 



/(7)d7 



(7/7o)3[l- (7/70)2^,-6] 



F2 = ^70/ [w^7o] , 

v.p. means that the integral is taken in the prin- 
cipal value sense, /(7) is the particle distribution 
function with the normalization / /(7)d7 = 1, 70 
the characteristic Lorentz-factor of the plasma par- 
ticles, Tc the characteristic radius of cyclotron res- 
onance defined as tj^f (rc) = uj^aO'^ /2. In equation 
(5) we have omitted the factors common to all the 
Stokes parameters, since they are irrelevant to the 
problem of polarization evolution considered. ^ The 
set of e quations f5) is a generalization of equations 
(17) in lPetroval ||200Cf) for the case of arbitrary ^. 
At the same time, we still do not include the terms 
responsible for rotational aberration explicitly and 
assume that they enter bx and by as factors of or- 
der unity (for more detail see Petrova 2006). The 
terms containing F\ describe polarization evolu- 
tion because of wave mode coupling, while those 
containing F2 correspond to cyclotron absorption. 
The initial conditions for the natural waves read: 
{I,q,u,v)o ~ (1, ±1,0,0), where the upper and 
lower signs refer to the ordinary and extraordinary 
waves, respectively. 

To proceed further, we choose a simplified dis- 
tribution function of the plasma particles in the 
form of a triangle: 



7 — 70 (1 ~ ") 

/(7)=^ 



70 



, " 70 
(1 + a) 



i i 
a 7o 



7 



70(1 - a) < 7 < 70 • 
7o < 7 < 70(1 + a) • 



^ Note, however, that the cyclotron absorption can 
markedly affect the total intensity of pulsar radiation. 
In the short-period pulsars, P ~ 0.1 s, the optical depth 
to th is process can exc eed uni t y (L vubarskii & Pctrov^ 
|1998?). As is shown in Petrova] 120021) . the effect of reso- 
nant absorption can account for the observed statistical 
features in the energetic characteristics of the short- 
period pulsars. 



Generally speaking, the numerical results on po- 
larization evolution are not very sensitive to the 
detailed shape of /(7) and are mainly affected by 
the average Lorentz-factor, 70, and the width of 
the distribution (the parameter a in equation (6)). 
The set of equations (5) along with equation (6) 
describe polarization evolution of the original nat- 
ural waves in the hot magentized weakly inhomo- 
geneous plasma of pulsars. 

Since the natural waves are completely polar- 
ized by definition, their polarization state is de- 
scribed by two independent quantities, e.g. PA, 
■(/i = 0.5 arctan(u/g), and the ellipticity, x = 
0.5 arctan " Figure 2 shows an example of 



(6) 



numerical tracings of tp and x of the original ordi- 
nary and extraordinary waves as they propagate in 
pulsar magnetosphere. One can see that polariza- 
tion evolution can indeed be significant and differs 
markedly for the two modes. 



3 STATISTICS OF THE FINAL 

POLARIZATION PARAMETERS 

3.1 Histograms of PA and ellipticity 

As is discussed above, polarization evolution is de- 
termined by the parameters /x and rj and is affected 
by the parameters of the particle distribution func- 
tion. For the rays observed at a fixed pulse longi- 
tude, the locations of the regions of wave mode 
coupling and cyclotron resonance may vary from 
pulse to pulse because of fluctuations in the num- 
ber density and the characteristic Lorentz-factor of 
the plasma particles. Besides that, the fluctuations 
in the plasma distribution may affect refraction of 
the waves, so that the rays of a certain orientation 
may follow somewhat different trajectories in the 
open field line tube. All this is believed to influence 
the final polarization states of the original natural 
waves. Note that although the modern theories of 
the pulsar pair creation cascade say nothing about 
the fluctuations in the resultant distribution of the 
secondary plasma, an idea of such fluctuations is 
supported by the random nature of the cascade 
process and is strongly suggested by the variabil- 
ity of the observed individual pulse proflles. 

In the upper panel of Fig. 2, the left plot shows 
the histogram of the final PA of the original natu- 
ral waves in case of randomly varying parameters 
fi and rj. One can see that the peaks are separated 
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by not exactly 90° , the scatter of PA values around 
the peaks is about 10° and the ordinary mode (that 
with positive PA) shows somewhat less scatter than 
the extraordinary one. The histogram of mode el- 
lipticities for the same fi and rj is shown in the 
right plot of the upper panel in Fig. 2. Note that 
the peaks are not symmetrical with respect to zero 
and the extraordinary mode exhibits a more pro- 
nounced scatter. 

The middle and bottom panels of Fig. 3 show 
the histograms of PA and ellipticity for the sum 
of the two polarization modes given random vari- 
ations in fi, rj and the mode intensity ratio. In the 
present consideration, we assume that initially only 
the ordinary mode is emitted, whereas the extraor- 
dinary one arises as a result of partial conversion 
of the ordi nary mode de ep inside the magneto- 
sphere (see iPetroval l200lh . This is supported by 
the recently discovered anticorrclation of the mode 
intensities (Ed wards fc S tappcrs 2004). Since the 
mode conversion is a propagation effect, its effi- 
ciency is also determined by the instantaneous dis- 
tribution of the plasma and is believed to fluctuate. 
Thus, the coefficient of conversion, r, is considered 
as a random quantity. Generally speaking, it may 
be correlated with n and r), but as the process of 
mode conversion takes place far from Vp and r^, 
the correlation is thought to be weak and is ne- 
glected throughout the paper. In the middle panel 
of Fig. 2, T is taken to be distributed over the whole 
interval from to 1, whereas in the bottom panel 
T G [0.6,1]. 

Although in the former case the probability to 
dominate is the same for the two types of original 
natural waves, the humps at the PA histograms 
are substantially distinct: the peak corresponding 
to the ordinary waves is much more pronounced. 
This is a consequence of cyclotron absorption, 
which suppresses the extraordinary constituent of 
the wave polarization stronger. The extraordinary- 
mode hump looks more smeared and is connected 
to the ordinary-mode one with a bridge. In the bot- 
tom panel, the humps change in dominance, but 
the ordinary mode is still present, despite com- 
plete dominance of the extraordinary waves just 
after conversion. This is again because of differ- 
ential action of cyclotron absorption on the two 
types of natural waves. The bridge between the 
humps looks more pronounced, and the histogram 



on the whole is similar to the observed ones (e.g. 
iMcKinn on 20031. 

In the middle histogram of the ellipticity, the 
hump corresponding to the extraordinary mode is 
barely resolved, whereas the ordinary mode peaks 
at X ~ 0. On the whole, this distribution looks 
like a unimodal one with a long tail. As can be 
seen in the bottom histogram of x, only positive 
values are met and the two humps are barely re- 
solved. This distribution can also be regarded as 
a unimodal one, in contrast to the corresponding 
histogram of PA. It is not our aim here to fit the 
concrete distributions observed, but Fig. 3 demon- 
strates the principal possibility of such fits within 
the framework of the propagation model of pulsar 
polarization. 



3.2 Two-dimensional scatter-plots 

Propagation origin of pulsar polarization implies 
a certain correlation between the mode elliptic- 
ity and PA: Both these quantities are determined 
by the instantaneous state of the plasma and vary 
from pulse to pulse because of fiuctuations in the 
plasma fiow. Now we are going to model the obser- 
vational signatures of this correlation for different 
sets of fj, and rj and compare our simulations with 
the observational results. 

At a fixed pulse longitude, the simulated 
individual-pulse polarization presents a sum of the 
two polarization modes with random intensity ra- 
tio. The resultant Stokes parameters are plotted 
in the Lambert's azim uthal equal area projection, 
as is done by observers jEdwards fc Stappersl20o3 : 
Edwards 20o3)- 

As each mode can dominate from 
time to time, at the Poincare sphere the resultant 
Stokes vectors tend to group around the roughly 
orthogonal locations corresponding to the observed 
modes. Let us choose the spherical coordinate sys- 
tem with the polar axis along the fiducial vec- 
tor of the normalized Stokes parameters, Sp = 
{qp,Up,Vp}, qp + Up + Vp = 1, which is the aver- 
age over the Stokes vectors belonging to one of the 
observed modes. The quantities p and A are the po- 
lar angle and azimuth of the normalized individual 
Stokes vectors (q,u,v) in this spherical system: 



p — aiccos(qqp + uUp -f vvp) . 



A = arctan 



uqp — qup 



{qqp + uup)vp — v{ql + u^) 



(7) 
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If one consider p and A as the polar coordinates 
(the radius and azimuth, respectively), we come to 
the Lambert's azimuthal equal area projection of 
the Poincare sphere. This projection is interrupted 
at the equator, and for the Stokes vectors lying in 
the southern hemisphere (i.e. for the second mode) 
the projection of the sphere with the polar axis 
along Sm = — Sp is conside red (for more detail see 
lEdwards fc Stapperjliooi) . 

Figure 4 shows the two-dimensional plots for 
the two original natural modes after their evolu- 
tion in the plasma in case of negligible cyclotron 
absorption, rj ~ 0.01, and moderate mode cou- 
pling, jj, ~ 0.1. Here Sp corresponds to the origi- 
nal ordinary mode. One can see that the process 
of mode coupling in the fluctuating plasma re- 
sults in the unambiguous relation of the polariza- 
tion characteristics, which has already been pro- 
posed as a bas is for diagnostics of pulsar plasma 
ilPetrovalliool) . At these plots, the points prefer 
certain azimuths, which is in qualitative agreement 
with the observational result of Edwards (2004) 
(see Fig. 5 there), though the scatter of the ob- 
servational points is enormously large. Note that 
this scatter cannot be reproduced by considering 
the sum of the two modes with random intensi- 
ties. Introducing a slight non-orthogonality of the 
modes because of cyclotron absorption (ry ~ 0.1) 
makes the plot more realistic (see Fig. 5 for the case 
of completely dominating original ordinary mode, 
r G [0, 0.3]), though the scatter of the points is still 
insufficient. 

Figure 6 shows the two-dimensional scatter 
plots in case of strong non-orthogonality of the 
modes (r; ~ 1). In the upper panel, the fiducial 
Stokes vector, Sp, corresponds to the state with 
the dominant ordinary mode, whereas the coeffi- 
cient of conversion is uniformly distributed in the 
interval [0.5, 1]. Note that although in this case 
the conversion of the ordinary waves into the ex- 
traordinary ones is efficient, the final polarization 
state of the sum of the modes can still be close to 
that of the ordinary ones because of less signifi- 
cant absorption of the ordinary waves and because 
of a marked non-orthogonality of the modes. The 
bottom panel of Fig. 6 corresponds to the case of 
a uniform distribution of r over the whole inter- 
val from to 1, with Sp directed along the ex- 
traordinary mode vector. The plots in both pan- 
els exhibit strong non-orthogonality of polarization 



states, with the secondary modes forming the char- 
acteristic arc-like features, which are in a quali- 
tative agreement wit h those discovered in obser - 
vations (cf. Fig. 2 in lEdwards fc Stapperj|2004^ . 
Note, however, that the scatter of the simulated 
points and the degree of non- orthogonality are sub- 
stantially less than those in the observational plots. 
This may be an indication of the role of some ad- 
ditional effects in the formation of the single-pulse 
polarization of pulsars. 



4 DISCUSSION AND CONCLUSIONS 

The hot magnetized weakly inhomogeneous plasma 
of pulsars may substantially afi'ect the radio wave 
propagation. The polarization states of the original 
natural waves may change markedly because of the 
wave mode coupling and cyclotron absorption. The 
former process turns the original linearly polarized 
waves into the elliptical ones, which are still purely 
orthogonal at the Poincare sphere. The mode cou- 
pling efficiency is determined by the location of the 
coupling region, rp, in the tube of open magnetic 
lines. The role of cyclotron absorption in the evolu- 
tion of wave polarization depends on how close the 
regions of mode coupling and cyclotron resonance 
are, being most prominent in case of their approx- 
imate coincidence, rj — Vp/rc ~ 1. Typically rc 
lies beyond rp, in which case cyclotron absorption 
leads to the non-orthogonality of the waves. Note 
that all the way in the magnetosphere the original 
natural waves remain completely polarized. 

Temporal fluctuations in the plasma flow are 
believed to underlie the fluctuations of the individ- 
ual pulse polarization. The plasma distribution in 
the open field line tube can be strongly variable 
because of the random character of the pair cre- 
ation process. Hence, for the rays observed at a 
fixed pulse longitude, both rp and rc can change 
from pulse to pulse due to variations of the num- 
ber density and characteristic Lorentz-factor of the 
plasma. Besides that, the variations of the plasma 
distribution may affect refraction of waves, so that 
the rays observed at a given pulse longitude may 
follow somewhat different trajectories in the mag- 
netosphere and have different tilts to the ambient 
magnetic field while passing through the regions of 
mode coupling and cyclotron resonance. 

The intensity ratio of the modes fluctuates 
as well, and it is also thought to result from the 



© 0000 RAS, MNRAS 000, 000-000 



8 S. A. Petrova 



plasma fluctuations in the magnetosphere. The 
point is that the extraordinary natural waves have 
vacuum dispersion and can hardly be generated 
directly by any conceivable emission process in 
the plasma. Therefore they are likely to originate 
as a result of partial conversion of the ordinary 
waves, which can take place in the regions of quasi- 
longitu dinal propaga tion deep inside the magneto- 
sphere JPetroval200j) . The idea of mode conversion 
is supported by the recently discovered anticorrcla- 
tion of the mode intensities (Edwards & StapDcrs 
l200dl . In the framework of this view, the mode in- 
tensity ratio is determined by the coefficient of con- 
version and changes from pulse to pulse because of 
fluctuations in the plasma flow. 

Thus, the propagation model of pulsar po- 
larization incorporates two superposed completely 
polarized modes with randomly varying polariza- 
tion states and intensities. It is important to note 
that because of cyclotron absorption the super- 
posed modes can become non-orthogonal, whereas 
the original natural waves are purely orthogonal by 
deflnition. 

In the present paper, the statistics of the indi- 
vidual pulse polarization have been simulated un- 
der the assumption that the parameters /i and rj 
as well as the coefficient of conversion, r, are the 
random quantities with unifirm distributions over 
some intervals. The resultant histograms of PA ex- 
hibit two humps markedly smeared and connected 
with a bridge. The peaks can be separated by not 
exactly 90°. The histograms of the resultant el- 
lipticity look like the unimodal ones because of a 
very small separation between the peaks of the two 
observational modes. Furthermore, it may happen 
that in the whole sample the ellipticity is purely of 
one sign, whereas the corresponding histogram of 
PA is bimodal. Although direct fits to the obser- 
vational data are beyond the scope of the present 
paper, our results confirm the ability of the prop- 
agation model to account for the main features of 
the observed histograms of PA and ellipticity. 

Propagation origin of pulsar polarization im- 
plies a certain correlation between the mode ellip- 
ticity and PA. Given that the contribution of cy- 
clotron absorption is negligible and the final polar- 
ization is determined solely by the mode coupling, 
there is a one-to-one correspondence between the 
mode ellipticity and PA. However, it is not proved 
by the observational data available. Although an 



evidence for the expected relation is indee d present 
in th e polarization of PSR B0818-13 (lEdwardsl 
120041) . the scatter of the observational points ap- 
pears dramatic, and it cannot be reproduced by 
taking into account the fluctuations of the mode in- 
tensity ratio. In case of moderately weak cyclotron 
absorption, a slight non-orthogonality of the fluctu- 
ating modes causes an additional scatter of the po- 
larization parameters and allows to reproduce the 
observational plots on a qualitative level, though 
the scatter is still insufficient quantitatively. In 
case of strong non- orthogonality of the modes, our 
model qualitatively reproduces the characteristic 
arc-like features prese nt in the observational plot s 
for PSR B0329-H54 jEdwards fc StappersI liooi) . 
though the magnitude of non-orthogonality and 
the total scatter of the simulated points are again 
less than in observations. 

Thus, on a qualitative level, the expected cor- 
relations of polarization parameters are compati- 
ble with the observational data. This is a strong 
argument in favour of the propagation model of 
pulsar polarization. Apparently, to achieve better 
quantitative agreement with the observational re- 
sults it is necessary to improve and further de- 
velop the model suggested. In the present consid- 
eration, we have ignored the net charge density 
in the magnetosphere and treated the process of 
cyclotron absorption in the limit of small pitch- 
angles of the particles. In reality, both these as- 
sumptions can be violated, so that our results on 
the individual-pulse polarization and its statistics 
can be somewhat modified. Besides that, pulsar 
emission can be considered as a su m of contri- 
butions from multiple sub sources (e.g. lGil fc Lvnel 
I1995I: iMelrose et ahllioo^ . in which case the pos- 
sibilities of modeling the resultant polarization are 
much wider. 

In conclusion, it should be pointed out that 
the model suggested allows a unique possibility 
of diagnostics of pulsar plasma by means of the 
individual-pulse polarization. 
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Figure 1. Geometry of the vectors fe and b in the coordinate system described in the text. So is the initial orientation 
of the magnetic axis. 
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Figure 2. Position angle (a) and ellipticity (b) of the original ordinary and extraordinary modes as functions of 
distance in pulsar magnetosphere; n = 0.3, »? = 1, a = 0.1. 
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Figure 3. Numerically simulated histograms of the final PA and ellipticity after polarization evolution in the fluctu- 
ating plasma. The parameters /i and rj are uniformly distributed over the intervals [0.2,0.4] and [0.9,1.1], respectively. 
The upper panel shows the histograms for the original ordinary aud extraordinary modes. The middle and bot- 
tom panels correspond to the sum of modes with random intensity ratio. The coefficient of conversion is uniformly 
distributed in the intervals [0,1] and [0.6,1], respectively. 
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Figure 4. Two-dimensional scatter plots of the final Stokes parameters of the original ordinary (left panel) and 
extraordinary (right panel) modes in the Lambert's aaimutlial equal area projection (for more detail see text); 
/U £ [0.1,0.5] and rj € [0.01,0.05]. The one-to-one correspondence of the polarization parameters is evident. 
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Figure 5. The Lambert's aaimuthal equal area projection of the final Stokes parameters of the sum of modes with 
the dominant ordinary mode; r 6 [0, 0.3]; the role of cyclotron absorption is not negligible, rj £ [0.1, 0.4]; /i 6 [0.1, 0.5]. 
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Figure 6. Two-dimensional scatter plots of the final Stokes parameters in case of strong non-orthogonality of the 
modes, rj G [0.5, 1.5]. In the upper and lower panels, r G [0.5, 1] and T € [0, 1], respectively, and Sp is centered on the 
ordinary- and extraordinary-wave positions, correspondingly. 



© 0000 RAS, MNRAS 000, 000-000 



